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In this article, we present an ongoing lattice study of the Higgs- Yukawa model, in the regime 
of strong- Yukawa coupling, using overlap fermions. We investigated the phase structure in this 
regime by computing the Higgs vacuum expectation value, and by exploring the finite-size scaling 
behaviour of the susceptibility corresponding to the magnetisation. Our preliminary results indi- 
cate the existence of a second-order phase transition when the Yukawa coupling becomes large 
enough, at which the Higgs vacuum expectation value vanishes and the susceptibility diverges. 



The XXIX International Symposium on Lattice Field Theory - Lattice 2011 
July 10-16, 2011 

Squaw Valley, Lake Tahoe, California 



♦Preprint number: DESY 11-219 
^ Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Strong-Yukawa coupling 



C.-J. David Lin 



1. Introduction 

In recent years, there have been interests in the possible existence of heavy extra-generation fermions 
(mass > 600 GeV) beyond the standard model (SM). Such heavy fermions are a consequence of 
strong-Yukawa couplings. Their presence in nature remains to be examined by experimental data 
that will be collected at the LHC. An important consequence of a 4th fermion generation is the 
substantial enhancement the amount of CP violation Large bare values of the Yukawa cou- 
phng may also lead to the formation of bound states which can replace the role of the Higgs boson 
in unitarising the WW scattering process |Q ^, Q]. Such a scenario is clearly of nonperturbative 
nature and motivates the use of lattice field theory as a first-principle and nonperturbative tool for 
this research avenue. Lattice investigations at small and moderate values of the bare Yukawa cou- 
phng showed ^ that the lower Higgs boson mass bound is strongly affected by the presence of 
a heavy 4th fermion generation when compared to results using a physical value of the top quark 
mass Still, in these simulations no signs of bound states were observed, as expected. 



However, lattice simulations have also revealed the existence of an interesting phase structure of 
the model at large values of the Yukawa coupling ^ ^ |T^, |n|, 12, 13, 14, 15]. These simulations. 



performed around 1990, were lacking an exact chiral symmetry on the lattice and the results of 
these works are therefore not easy to interpret and to connect to the SM. Recently, exact lattice 



chiral symmetry [16] were established and, in fact, lattice simulations employing this lattice chiral 



symmetry confirmed the phase structure at strong bare Yukawa coupling [ ]17| , |1 8[ ] as found in the 
earlier studies. 

In this work we further explore the phase structure of a lattice chirally invariant Higgs-Yukawa 
model at large values of the Yukawa coupling. Our main aim is to start a systematic investigation, 
whether the phase transitions between the symmetric phase with vanishing vacuum expectation 
value (VEV) v = and the broken phase with v > are governed by critical exponents that dif- 
fer from the (Gaussian) one of the SM. This is a highly non-trivial question. In particular, it 
can be shown that the lattice Higgs-Yukawa model in the limit of infinite bare Yukawa coupling 
reduces to a pure scalar non- linear a-model lH 0, |l^] with again Gaussian critical expo- 

nents. Hence, the most interesting and important question is, whether at large but finite value of 
the Yukawa-coupling a new and non-trivial universality class emerges. We emphasise that a quan- 
titative determination of critical exponents of the phase transitions in the strong Yukawa coupling 
region using chiral invariant lattice fermions was never attempted before. The above potential of 
nonperturbative physics at large Yukawa couphngs motivates clearly such an investigation. 



2. Simulation details 

The discretisation of the scalar- field theory leads to the action (with lattice spacing a set to 1) 



5, = -l9>.%+I 



(2.1) 
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where a labels the four components of the scalar fields, d is the number of the space-time dimen- 
sions, mo is the bare mass and Ao is the bare quartic self-coupling. For practical lattice simulations, 
it is convenient to perform the change of variables 



(P = V2k^, mt) = , ^0 = 7:2, <^2.2) 

where K is the hopping parameter. This renders the lattice scalar field theory to the Ising form 



5^ = -2^£o,%+£ o."+A(o."-i: 



(2.3) 



which is more suitable for exploring the phase structure. 

In this work, we use the overlap operator ^(°^) in the lattice fermion action 

Sf = where ^ = +P+ O'^ diag(3;,/,3;fo/) P++P- diag(3;,/,3;fa/) O P , (2.4) 

with 

^^t)-^ -(:»-::')■ 

The chiral projectors are defined as 



where p is the radius of the circle of eigenvalues in the complex plane of the free overlap operator. 
We also set yf = = y, to ensure that the fermion determinants are positive definite. 
Our simulations have been performed at two k values, K = 0.00 and 0.06, with the bare Yukawa 
coupling y in the range between 14 and 25. The bare scalar quartic coupling X is fixed to infinity, 
which results in the largest possible Higgs mass |^, |6[ 0]. We already accumulated data with reason- 
able statistics for the volumes 8^ x 16, 12^ x 24 and 16^ x 32. In generating dynamical scalar field 
configurations, we use the polynomial HMC (pHMC) algorithm [|T9|], treating the weight factor as 



an observable [20]. We have found that high degrees of polynomials {r^ 180 for 16^ x 32 lattices) 
are necessary in order to obtain reasonable statistical accuracy for the weight factor. For each set 
of the bare couplings, we carry out 1000 pHMC trajectories to precondition the fermion matrix 
and to thermalise the simulation. Our measurements are then performed on ~ 2000 thermalised 
trajectories. 

3. The scalar vacuum expectation value 

We first measure the scalar VEV to probe the phase structure. Our simulations have been performed 
without external sources, therefore a naive computation of the VEV would always lead to vanishing 
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Figure 1: Bare scalar VEV at two K values with X = °°. 



results in finite volume. In this work, we follow the procedure in Refs. |21, 22, 23] to project the 
scalar fields on the direction of magnetisation 



m : 




(V4 is the 4— dimensional volume), 



(3.1) 



then compute the VEV. This "projected" VEV coincides with the scalar VEV in the infinite-volume 
Umit, and its introduction in finite volume is equivalent to coupling the scalar fields to external 



sources |23] 



The results for the bare VEV, computed from the above projection procedure, are shown in Fig. [T]^ 
It is clear that there is a phase transition when y becomes large, at which the system enters a sym- 
metric phase. These plots also indicate that the value ofy = jcrit at which the phase transition occurs 
grows with k:. This agrees with the qualitative predictions from the strong-coupling expansion [^], 
the large-Nf expansion and an exploratory numerical study using overlap fermions [18]. 



4. Finite-size scaling of the magnetisation susceptibility 

In order to determine the order of the phase transition observed in the last section, we investigate the 
finite-size scaling behaviour of the susceptibility corresponding to the magnetisation. It is defined 
as 

X=V4{{m^)-{m){m)), (4.1) 



where V4 and m are defined in Eq. (3.1 ). This quantity diverges at the critical points in the infinite- 
volume limit. Finite-size effects result in the crossover from this bulk behaviour to the finite- 
volume scaling behaviour in lattice calculations. In the vicinity of a would-be second-order phase 



'We are currently computing the Goldstone wavefunction renormalisation to obtain the renormalised Higgs 
VEV lEl. 
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Susceptibility, K=0.06 
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Figure 2: Susceptibility at two K values at various volumes with X =oo. 



transition, the solution of the renormalisation group equation (RGE) in finite volume predicts [ ]26| ] 



r/v 



g(fLy^), with t = {y/{y,n, -A^/L'l) - 1), 



(4.2) 



where g is a universal function, is the spatial volume, jcrit is the critical Yukawa coupling at 
which the phase transition occurs in the infinite-volume limit, A4 is a phenomenological parame- 
ter, V and 7 are the universal critical exponents (anomalous dimensions), and b is the shift expo- 



nent ||27|]. 

For each of the two K values in our simulations, we perform a simultaneous fit of the data for the 



susceptibility at all volumes, to the partly-empirical formula [28] 



y-jcni-M/L, 



-y/2 



(4.3) 



where A 1,2,3,4 are unknown phenomenological coefficients. They are determined, together with v, 
7> Jcrit and b, from the fits. For our best procedure, the fit ranges of y are (14.5, 19.5) for = 0, 
and (14,22) for K = 0.06. The extracted 3'crit> 7> V and b are presented in Table [T[ For comparison 
with the scaling behaviour in the regime where the Yukawa couplings are weak and perturbative, 
we also list the predictions from the mean-field calculation in the 0(4) scalar model. We use these 
fit results to obtain the susceptibility as a function of y according to Eq. ( [4.3[ ). This is plotted with 
our data points in Fig. We also use the same fit results to construct x^-s ^^^^and tL}J^ , to examine 
the finite-size scaling behaviour of Eq. (12). The outcome of this test is shown in Fig. ||. 



Our data, as presented in Figs. § and ^, establish evidence for the existence of a second-order phase 
transition in the strong-Yukawa coupling regime. To further investigate the nature of this strong- 
Yukawa symmetric phase, we have to perform more detailed studies on the critical exponents and 
the spectrum. As demonstrated by the results collected in Table 7 is almost consistent with the 
corresponding mean-field prediction for the 0(4) scalar model, while v exhibits a more significant 
deviation from it. By varying the fit ranges in y in a reasonable interval for the above finite-size 
scaling analysis, the shifts in jcrit and the critical exponent 7 are not statistically distinguishable. 
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On the other hand, V can change by ~ 1 1 % and become almost consistent with the value in the 
0(4) scalar model. We have also tried incorporating logarithmic volume effects in Eq. (43) by 
replacing the coefficient A4 with A4[l +c Log(Lv)] (c is an unknown parameter). This procedure 
leads to no statistically significant variation of the critical exponents, V and 7. 



5. Summary and outlook 

In this article, we present an ongoing numerical investigation of the phase structure of the strong- 
Yukawa model on the lattice. Using overlap fermions which respect exact lattice chiral symmetry, 
and performing simulations at various volumes, enable us to explore the details of the phase struc- 
ture. From our computation of the scalar VEV and the study of the finite-size scaling behaviour 
of the magnetisation susceptibility, we obtain strong evidence that there exists a symmetric phase 
in the strong- Yukawa coupling regime, and that the transition between this phase and the broken 
phase is of second-order nature. At the moment, we cannot determine if the critical exponents for 
this phase transition are different from those in the weak- Yukawa regime. 

We are currently generating large lattices (24^ x 48) which will allow us to have more precise 
extraction of the critical exponents. These lattices will also enable us to control the infinite-volume 
extrapolations in our future calculation for spectral quantities. 





K = 0.00 


K = 0.06 


0(4) scalar model 


Vcrit 


16.57 ±0.06 


18.11 ±0.06 


N/A 


7 


1.02 ±0.02 


1.08 ±0.01 


1 


V 


0.57 ±0.03 


0.66 ±0.02 


0.5 


b 


2.05 ±0.20 


2.04 ±0.20 


N/A 



Table 1: The critical Yukawa coupling jcrit, the critical exponents 7, v, and the shift exponent b determined 
from the best fits to the finite-volume scaling function. The errors are statistical only. Predictions from the 
mean-field calculation in the (9(4) scalar model are also listed. 
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